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Summary.  Annual  mean  sunspot  numbers  R(t)  since  1700  show  evidence  of  a  nonlinear  effect,  first 
evidenced  by  the  detection  of  third  harmonic  in  R^t),  the  alternating  representation  of  the  magnetic 
(22-year)  cycle  of  solar  activity.  The  form  of  the  nonlinearity  proves  to  be  a  three-halves  law  R(t)  = 


(1 


mRun(t)) /83jf^where  Ah*(0.  also  an  alternating  quantity,  is  a  presumed  underlying  or  '‘linearised’’^ 
sunspot  number.  The  nonlinearity  is  of  such  a  nature  as  to  cause  strong  semicydes  to  be  sharper  than 


sinusoidal  and  to  produce  the  inflexion  in  noted  at  sunspot  minimum.  The  difference  R(t)  —  |/2j,n(t)| 
is  sufficiently  like  third  harmonic,  for  semicycles  of  average  strength,  tO'fexplain  the  band  around  22/3 
years  which  is  noticeable  in  the  spectrum  of  R±(t).  However,  third  harmonic  alone  is  not  sufficient  to 
account  for  the  observed  dependence  of  eemkyde  shape  on  amplitude,  whereas  the  three-halves  law  accounts 
economically  for  a  range  of  effects^  A  search  for  a  physical  explanation  of  a  three-halves  law  reveals  that 
just  such  a  law  results  because  large  sunspot  groups,  such  as  occur  around  strong  sunspot  maxima,  enter 
into  the  sunspot  number,  as  conventionally  defined,  over  more  days  than  small  groups,  simply  because  large 
groups  last  longer. 

Semicycle  asymmetry,  which  cannot  result  from  a  simple  nonlinear  law,  is  here  ascribed  to  magnetic 
buoyancy  acting  preferentially  on  the  antinod al  layers  of  a  travelling  wave.  Profiles  for  semicydes  of  different 
strengths  have  been  constructed  on  the  assumption  that  the  underlying  influence  is  sinusoidal.  Each  sinusoid 
is  distorted  by  the  three-halves  law,  and  then  made  u asymmetrical  by  applying  a  buoyancy  theory  for 
magnetised  plasma  rising  against  viscous  drag.  The  majority  of  past  semicydes,  including  those  of  the 
18th  century,  can  be  matched  quite  well  by  the  artificial  profiles,  a  conclusion  that  supports  the  idea  of  an 
underlying  influence  that  is  sinusoidal,  and  the  hypothesis  that  sunspots  result  from  an  upward  traveling 
wave  from  a  submerged  22-year  oscillator. 
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Introduction 

The  shape  of  the  sunspot  cycle  has  been  of  importance  for  many  years  in  connection  with  ionospheric 
physics,  magnetospheric  physics,  and  interplanetary  phenomena,  including  solar  proton  events,  partly  be¬ 
cause  of  the  impact  on  such  diverse  practical  matters  as  the  scheduling  of  high-frequency  radio  communica¬ 
tion,  over-the-horison  radar,  magnetic  mines,  ephemerides  of  low-level  satellites,  and  the  safety  of  astronauts. 

The  eleven-year  recurrence  tendency  has  been  variously  attributed  to  an  oscillator  subject  to  distur¬ 
bances  or  alternatively  to  a  resonator  excited  by  random  impulses;  one  sees  that  these  proposals  range  from 
an  active  source  at  the  one  extreme  to  a  passive  structure  at  the  other.  Among  mathematical  models  pro¬ 
posed  in  the  past,  some  have  been  based  on  these  two  not  very  definite  but  nonetheless  physical  approaches, 
while  other  models  have  been  frankly  empirical  with  no  physical  basis  at  all.  An  understanding  of  the  shape 
of  the  solar  cycle  would  be  helpful  for  prediction  in  practical  circumstances. 


Features  to  explain 

The  traditional  solar  cycle  begins  at  sunspot  minimum,  lasts  about  11  years,  and  ends  at  the  following 
minimum,  but  since  magnetic  polarities  reverse  at  sunspot  minimum,  an  ambiguity  has  developed  in  the 
meaning  of  the  word  cycle.  For  example,  one  may  read  that  the  solar  ‘‘cycle’’  is  really  22  years.  For  this 
reason  the  terms  “semicycle”  and  “magnetic  cycle”  will  be  used  to  distinguish  the  11  ±  1.5b- year  episodes 
from  the  22.2  ±  2.08-year  ones. 

The  course  of  an  *1 1-year"  semicycle,  as  represented  by  annual  means,  has  a  distinctive  nonsinusoidal 
shape  that  varies  from  occasion  to  occasion,  apparently  in  dependence  upon  at  least  two  parameters,  the 
duration  T  and  the  maximum  R^ot  reached  by  the  annual  mean  sunspot  number  R.  Starting  from  sero, 
or  a  small  value,  at  (  =  0,  there  is  an  initial  fillet  where  tPR/dt *  >  0,  a  convex  segment  peaking  at  t  =  tm, 
followed  by  a  more  or  less  gentle  decline  toward  sero  at  t  =  T  that  is  often  separated  from  the  convex 
segment  by  a  noticeable  break  in  downward  slope.  Usually  the  maximum  occurs  earlier  than  at  t  =  J7\ 
the  more  so  as  A*,*,  is  larger;  roughly  speaking,  the  product  Rm„tm  tends  to  be  constant.  In  contrast, 
the  declining  phase  tends  to  be  longer  as  /2mU  is  larger,  the  ratio  Rmao/(T  -  tm)  tending  to  be  constant. 
The  unsymmetrical  aspects  are  leas  pronounced  when  the  semicycle  is  a  weak  one.  In  Fig.  1,  on  the  left, 
these  features  are  exhibited,  for  the  strong  semicycle,  by  the  average  of  the  semicydes  beginning  1775,  1784, 
1833,  1843,  1867,  1933,  1944,  1954,  1976,  all  of  which  had  Rmao  >  114.  The  shape  of  a  weak  semicycle  is 
represented  by  the  average  for  1755,  1798,  1810,  1823,  1878,  1889,  1901,  1923,  where  R^  <  86. 


•  Fig.  1.  Features  of  the  shape  of  a  strong  semicycle  (left)  indude  a  leading  fillet,  a  convex  segment 


peaking  early,  and  a  relatively  long  decline,  sometimes  with  a  noticeable  break  in  slope.  A  weak  semicycle 
(right)  tads  to  be  much  more  symmetrical. 


Fourier  analysis,  derectification  and  linearisation 

Fourier  analysis  of  the  annual  mean  sunspot  number  series  R(t)  has  been  frequently  carried  out  but 
it  has  always  been  difficult  to  interpret  the  results.  If  one  agrees  that  R(t )  =  |A±(f)|,  where  R±(t)  is  an 
alternating  series  that  reverses  sign  in  successive  semicycles  (Brace well,  1953)  and  whose  period  is  the  full 
magnetic  cycle,  then  it  is  easy  to  see  where  much  unsignificant  spectral  detail  has  originated:  R(t)  is  related  to 
A±(t)  in  the  same  way  that  the  output  voltage  of  a  full- wave  rectifier  is  related  to  its  input;  thus  rectification 
of  the  alternating  series  R±(t)  results  in  the  appearance  in  R(t)  of  sum  and  difference  frequencies  without 
physical  relevance.  Fourier  analysis  of  by  contrast,  gives  a  much  simpler  spectrum  (Cole,  1973).  The 
factitious  band  of  spectral  lines  around  11. 1-1  c/a,  together  with  the  harmonics  of  this  band,  practically 
disappears,  while  a  cleaner  22.2-1  c/a  band  is  sea.  Not  only  does  the  naive  spectrum  of  R(t)  contain 
artifacts  of  rectification,  but  other  significant  spectral  features,  such  as  a  third  harmonic  (period  ~  22/3 
years)  which  appears  in  Cole’s  spectrum,  are  suppressed.  This  particular  feature  has  recently  assumed 
considerable  physical  importance. 

Such  a  third  harmonic  was  independently  noticed  in  the  course  of  a  attempt  to  determine  an  envelope 
for  A±(i)  (Bracewell,  1985),  where  it  was  found  that  removal  of  the  harmonic  by  filtering  was  an  essential 
prerequisite  to  access  of  the  Hilbert  transform  envelope  procedure.  Since  then  a  similar  third-harmonic 
phenomenon  was  noticed  in  the  application  of  Hilbert  transform  analysis  to  varve  data  (Bracewell  and 
Williams,  1986).  For  the  purpose  of  those  two  papers  it  was  sufficient  to  reduce  the  unwanted  harmonic 
component  to  a  tolerable  level  by  filtering;  but  it  was  observed  that  the  waveform  being  suppressed  was  not 
strictly  a  monochromatic  third  harmonic  of  the  22.2-year  fundamental.  Instead,  the  sero- crossings  of  the 
unwanted  perturbation  were  so  related  to  each  “11- year"  semicycle  that  the  perturbation  was  aacribable  to 
a  nonlinear  effect  operating  year  by  year  on  the  oscillating  quantities,  annual  sunspot  number  ad  varve 
thickness  respectively.  Strong  semicycle*  were  accompanied  by  a  perturbation  of  amplitude  greater  than  in 
direct  proportion  to  the  strength  of  the  semicycle,  while  the  sero- to- sero  spacing  of  the  perturbation  was  not 
fixed,  as  would  be  true  of  a  harmonic,  but  changed  in  accordance  with  the  duration  of  that  semicycle.  Thus 
the  phenomenon  which  originally  revealed  itself  through  the  appearance  of  a  7-year  peak  in  the  computed 
spectrum  received  a  compact  description.  The  spectral  description,  on  the  other  had,  while  mathematically 
equivalent,  proved  to  be  dependent  upon  the  duration  ud  epoch  of  the  time  segment  chosen  for  spectral 
aalysia. 

Since  both  the  sunspot  series  ud  the  varve  series  exhibit  this  phenomenon  of  cycle  shape  distortion,  a 
source  in  the  an  is  indicated. 

Subsequently  it  was  found  that  the  sunspot  number  curve  could  be  retrodicted  by  the  introduction  of  a 
quadratically  nonlinear  function 

R±  =  (Ran  +  0Rfin)  »gn  Riin 

in  which  the  coefficient  was  adjusted  empirically  at  0  =  0.0035  (Bracewell,  1986).  If  the  linearised  sunspot 
number  Ru»  is  taken  to  vary  quaaisinusoidally  then  the  enhancement  at  large  R  will  cause  the  sharpened 
profile  familiar  in  strong  semicycles  of  observed  sunspot  number.  Linear  Fourier  synthesis  with  inclusion  of 
third  harmonic  is  unable  to  mimic  the  observed  semicycle  shapes  as  well  as  the  nonlinear  formula  allows. 

The  motivation  for  nonlinearity  correction  was  that  Ran  might  be  in  some  sense  physically  more  fun¬ 
damental  thu  the  anual  sunspot  number  JfZ(f),  a  presumption  that  is  supported  to  the  extent  that  some 
features  of  R(t)  are  accounted  for  ocTthe  spectrum  becomes  simpler. 

The  form  of  the  nonlinear  law 

Instead  of  aaaming  the  nonlinear  relation  represented  by  a  srbitary  quadratic  law  with  the  one  dispos¬ 
able  parameter  0  it  would  be  nice  to  determine  the  law  explicitly  point  by  point.  Suppose  that  measurements 
y(t)  were  available  of  a  nonlinear  function  H  of  a  quaaisinusoidally  varying  function  *(t)  =  A  sin( wt  +  o),  i.e. 
V<<)  =  7f(x(t)].  The  situation  could  be  repreauted  as  in  Fig.  2.  Now  the  function  “H  could  be  determined 
by  plotting  y{t)  versus  *(1)  with  1  as  a  parameter,  if  r(f)  were  known. 

•  Fig.  2.  A  nonlinear  characteristic  (top  left)  relates  a  stimulus  *(t)  (lower  left)  to  a  response  p(f) 
(upper  right).  In  terms  of  the  corresponding  points  (x*,,tm)  ad  (p*,,!*,)  the  reduction  factor  K  is  Sw/nu- 

la  this  preerat  problem  *(t)  is  not  known,  we  only  know  y(t),  ad  it  would  seem  hopeless  to  ask  for 
H.  But  we  do  know  that  A  ud  a  are  slowly  varying  functions  of  time,  a  property  that  provides  a  toehold. 
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A  remarkable  iterative  approach  haa  been  found  that  makes  use  of  the  quadratic  approximation  to  get  an 
immediate  improvement  in  the  law.  In  Fig.  2  suppose  that  y(f)  is  drawn  in  but  that  z(f)  and  the  desired 
law  H  are,  to  begin  with,  missing.  We  can  start  by  marking  the  seto  crossings  of  *(<)  on  the  down-going 
time  axis,  and  we  can  also  mark  the  direction  of  crossing.  This  is  because  we  are  looking  for  a  monotonic 
law  that  passes  through  the  origin.  The  slopes  z'(i)  at  z(t)  =  0  are,  however,  not  known,  contrary  to  the 
assumption  of  limiting  linearity,  built  into  the  quadratic  law,  according  to  which  R  is  proportional  to  Ran 
at  small  values.  If  the  law  were  linear  at  small  sunspot  numbers,  as  the  quadratic  expression  implies,  then 
the  slope  of  z( I)  at  a  sero  crossing  would  be  ±1.  Now,  however,  we  would  like  to  avoid  this  prejudgment 
and  determine  the  actual  law  from  the  data.  The  reasoning  proceeds  as  follows. 

In  the  time  interval  [-4/v,  (*— 4)/u>]  between  adjacent  sero  croesings  let  p(t)  have  an  observed  maximum 
Pm.  The  artificial  half-period  sinusoid  Pm  sm(w<  +  4)  can  be  constructed  and  compared  in  shape  with  p(t)  in 
that  time  interval  to  give  some  information  about  the  nonlinearity  in  the  range  0  to  If  such  information 

from  all  the  semicycles  could  be  combined,  knowledge  of  the  nonlinear  characteristic  could  be  sharpened.  To 
make  results  from  all  the  semicycles  superposable,  imagine  a  graph  of  p(t)  against  Pm  sin(u/t  +  4)  in  which 
the  abscissae  for  any  one  semicycle  are  reduced  by  the  constant  factor  K  =  [-1  +  (1  +  <fipm  )^*]/2/Jy^  for 
that  semicycle.  The  reduction  factor  comes  from  solving  pm  =  z  +  fiz*  and  dividing  by  Pm-  Superposing 
points  from  the  286  available  years  since  1700  would  confirm  the  quadratic  expression,  if  it  were  correct, 
because  all  the  points  would  fall  on  the  one  curve.  Otherwise  a  systematic  departure  would  be  seen  which 
would  be  the  basis  for  a  corrected  law. 

In  the  simple  form  described  this  approach  fails  on  the  sunspot  data  because  of  the  relation  between 
rise  time  and  maximum  discovered  by  Waldmeier,  as  a  concomitant  of  which  the  sunspot  maximum  does 
not  necessarily  fall  midway  between  sunspot  minima.  However,  given  the  years  (i  and  (j  of  two  adjacent 
sero  crossings  and  fm  the  year  of  maximum,  one  can  set  up  a  modulated  time-varying  phase 


d(0  =  a,(<  -  <i)  +  <*j(f  -  li)*, 


such  that 


One  can  show  that  the  coefficients  ate  given  by 


t  =  U 


o»  =  »(<«.  -  $(<i  +  h)]/(h  ~  -  t«)(U  -  <i) 

Oi  =  w/2(Um  -  <l)  -  OTj(U  -  ft). 


The  points  in  Fig.  3  show  the  result  of  plotting  R  ±  (t)  against  [-1  +  (1  +  4ffJfniu)£] 
nn4(i)/2ff  using  applicable  values  of  Rjim.  and  b(t)  for  each  semicycle.  Since  the  diagram  will  be  taken  as 
the  basis  for  a  new  law  defining  R^n,  the  abscissa  z  is  labeled 

•  Fig.  3.  The  noolinear  relation  between  R  ±  (()  and  Am»(<)  put  in  evidence  by  yearly  points  from 
1700  to  1985.  The  curve  is  R  ±  (I)  =  100{|  R,im(t)  |  /83)3'Jsgn 


Inflexion  effect 

The  rather  definite  functional  dependence  put  in  evidence  by  this  first  step  differs  only  a  little  in  general 
magnitude  from  the  quadratic  approximation  (which  is  not  shown),  but  has  a  noticeable  qualitative  difference 
at  low  values  of  sunspot  number.  One  notices  a  distinct  inflexion  rather  than  a  linear  passage  through  the 
origin.  Following  sunspot  minimum  the  rate  of  growth  of  sunspot  number  accelerates  for  a  year  or  two,  a 
feature  that  is  not  normally  thought  of  as  related  to  the  decline  of  the  sunspots  of  the  preceding  semicycle. 
In  the  alternating  view  afforded  by  R±(t),  the  phenomenon  appears  as  an  inflexion  in  what  would  otherwise 
be  a  more  or  leas  linear  crossing  through  sero.  The  rate  of  onset  of  a  new  semicycle  should  then  correlate 
with  the  rate  of  decline  of  the  old,  a  correlation  for  which  there  was  indeed  evidence  (Bracewell,  1953)  but 
which  the  inflexion  effect  tended  to  obscure. 

Curve  fitting 

lb  fit  a  curve  through  the  points  of  Fig.  3  is  a  classical  problem.  In  this  case  a  single  sided  logarithmic 
plot  (Fig.  4)  tells  us  that  a  simple  power  Taw  will  do  very  well.  With  this  knowledge  one  could  return  to 
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nonlogarithmie  coordinates  and  determine  parameters  by  least  squares.  However,  it  is  not  likely  that  one 
could  do  better  than  the  straight  line  fit  shown  on  Fig.  4  which  leads  to  the  revised  law 

R  =  100(I*«»I  /83)s/*  sgn  R,in. 

The  curve  drawn  on  Fig.  3  has  been  constructed  from  this  equation  and  clearly  is  about  as  good  a  fit  as 
one  could  expect,  given  the  spread  of  points,  although  in  principle  further  iteration  could  be  carried  out. 
The  straight  line  logy  =  1.5  logs  +  0.88  is  rather  sharply  defined;  if  the  three-halves  exponent  is  changed  by 
more  that  ±  0.1  an  acceptable  fit  cannot  be  found.  Consequently  for  practical  purposes  of  data  analysis  we 
shall  adopt  the  exponent  3/2. 

•  Fig.  4.  The  data  of  the  previous  figure  combined  in  one  quadrant  by  plotting  log  |y|  against  log  |ar|. 
The  straight  line  has  a  slope  3/2. 

The  method  of  solution  of  this  technical  problem  has  proved  to  be  most  satisfactory  and  could  be 
productive  in  applications  to  other  astronomical  and  geophysical  time  series. 

Effects  of  linearisation 

To  see  the  effect  on  sunspot  semicycle  shapes  that  is  produced  by  the  three-halves  law  consider  Fig.  5  in 
which  four  artificial  sinusoids  of  representative  strengths  (middle  row)  are  subjected  to  distortion  to  generate 
simulated  sunspot  semicycles  (top  row).  As  the  semicycles  become  stronger  they  develop  leptokurtosis,  or  a 
departure  in  the  direction  of  greater  sharpness  from  the  “normal”  condition,  taken  as  sinusoidal.  In  addition, 
the  inflexion  effect  in  H  at  low  values  strongly  influences  the  shape  of  the  weakest  semicycle  to  produce  a 
profile  reminiscent  of  the  records  obtained  around  the  years  1803,  1818,  1883  and  1907.  These  shapes  could 
be  described  as  platykurtic  in  the  dictionary  sense  of  the  word.  The  thing  that  distinguishes  the  three-halves 
law  from  the  quadratic  law  is  that  the  latter  produces  only  leptokurtosis  whereas  the  three-halves  law  gives 
rise  to  platykurtosis  also,  exactly  as  observed  in  the  annual  mean  sunspot  record. 

•  Fig.  5.  Artificial  “11-year"  semicycle  shapes  versus  time  (top)  associated  with  fixed-duration  semis- 
inuoids  (center)  of  four  different  strengths.  The  artificial  shapes  are  derived  from  a  three-halves  nonlinear  law 
without  introduction  of  asymmetry.  The  difference  curves  (below)  resemble  third  harmonic,  in  the  middle 
range  of  strength. 

The  difference  curves  (bottom  row)  show  that  the  perturbation  is  rather  like  third  harmonic  for  semicycle 
strengths  around  AmU  =  100  and  one  can  understand  how  Fourier  analysis  of  R  ±  (t)  could  show  a  band  of 
spectral  content  around  a  period  of  22.2/3  years.  However,  the  shapes  seen  for  very  weak  and  very  strong 
semicycles  show  that  linear  analysis  into  harmonics  does  not  do  justice  to  the  phenomenon  at  hand. 

Waldrueier  effect 

For  the  purposes  of  Fig.  5  symmetrical  sine  waves  were  used,  and  so  the  asymmetry  that  sets  in  with 
strong  semicycles  was  neglected.  Waldmeier  expressed  the  asymmetry  quantitatively  in  the  form  of  a  relation 
between  the  rise  time  and  the  maximum  sunspot  number  reached.  For  the  present  purpose  it  is  convenient 
to  express  the  asymmetry  in  terms  of  how  early  the  peak  arrives  with  reference  to  the  midpoint  between  zero 
crossings.  In  previous  work  (Bracewell,  1988)  this  lead  time  was  taken  as  1.5(/2m.,/150)1  years,  a  simple 
empirical  rule  that  is  quantitatively  compatible  with  Waldmeier’s  results.  However,  this  rule  makes  the  lead 
time  accelerate  with  increasing  Rm—%  which  is  unreasonable  because  the  lead  time  must  have  an  upper  limit 
(of  half  a  semicycle).  The  deficiency  becomes  apparent  when  sets  of  semicycle  shapes  are  calculated;  the 
shapes,  for  large  Rm—,  do  not  look  right.  Therefore,  a  dependence  based  on  a  theory  of  the  asymmetry  has 
been  attempted,  to  replace  the  empirical  relation.  The  following  derivation  is  based  on  the  idea  (Bracewell, 
1985)  that  asymmetry  is  due  to  magnetic  buoyancy.  The  only  other  explanation  previously  offered  for 
Waldmeier  effect  is  in  terms  of  the  eruption  hypothesis,  according  to  which  an  initial  impetus  carries  the 
sunspot  number  to  a  greater  or  lesser  height  from  which  it  recovers  with  a  time  constant  that  is  the  same 
for  all  semicycles.  This  is  more  descriptive  than  explanatory  and  lacks  a  physical  basis. 

Buoyancy  theory  for  semicycle  shapes 

Consider  an  oscillatory  source  of  horizontal  azimuthal  magnetic  field  H+  at  a  radius  ri  inside  the  sun, 
which  launches  waves  upward  with  period  22  years  and  a  phase  velocity  v.  At  a  greater  radius  r,  nodes  of 
arrive  with  a  delay  appropriate  to  the  velocity  v,  but  the  antinodes  arrive  early  because  of  magnetic  buoyancy. 


t 


4 


According  to  the  concept  of  magnetic  buoyancy,  introduced  by  E.N.  Parker,  solar  plasma  permeated  by  a 
magnetic  flux  density  B  experiences  a  vertical  buoyancy  force  proportional  to  B7,  because  the  magnetic 
pressure  relieves  the  particle  number  density  needed  for  quasihydroetatic  equilibrium.  We  now  apply  this 
concept  to  very  long  period  wave  motion,  keeping  in  mind  that  a  central  wave  source  is  onlv  a  working 
hypothesis  (Brace well,  19871,  and  that  general  opinion  tends  to  favour  the  relaxation  oscillation  originally 
introduced  by  Babcock  (19611  and  subsequent  developments  invoking  sonal  vortices  in  the  convective  region  •iJ. 
$nd  surface  dynamo  action  (Wilson,  1987). 

An  initially  sinusoidal  wave-field  variation  Bo  tin(fcr-wf)  in  the  vicinity  of  r  =  ri  should  suffer  distortion 
of  its  sinusoidal  form  as  it  slowly  propagates  upward.  The  general  effect  will  be  that  the  crest  will  rise  faster 
relative  to  the  preceding  and  following  nodes,  giving  the  waveform  a  forward-leaning  appearance  as  required. 

Given  an  initial  semicycle  profile  shape  Bi(r)  near  r  =  1*1  at  t  =  tj,  we  now  ask  for  the  distorted  shape 
near  the  greater  height  rj  at  a  later  time  t  =  ta. 

The  following  discussion,  which  is  offered  in  lieu  of  *  fall  theory,  assumes  that  the  22-year  wave  period  is 
so  long  that  buoyancy  forces  can  assert  a  quasiequilibrium  state  within  any  one  half-wavelength  layer.  From 
conservation  of  magnetic  flux  we  expect  that  B\(r)dr  =  (rj/ri)B2(r')  or 


The  very  long  time  scale  means  that  viscous  drag  dominates  inertial  forces  in  the  radial  direction.  The 
flux  element  B(r)dr  might  then  rise  with  an  excess  velocity  proportional  to  pressure  [Bi(r)]J.  Thus  r'  = 
r  +  ra  -  1*1  +  e[Bi(r))  ,  where  e  is  a  constant,  and  dr'/dr  =14-  2cBi(r)BJ(r).  So 


Bj(r)  = 


n  B>(r) 

r j  1  4-  2cBi(r)B[(r) 


Taking  Bi(rt)  =  A  sin[k(r  -  )] ,  ri  <  r  <  ^  +  JA,  we  find  that  an  initially  sinusoidal  shape  becomes 


B,(r)  = 


n _ Asin[k(r-  ra)] _ 

rj  1  4-  2cA  sin[k(r  -  r3)tA  cos [k(r  -  ra) 


rx  Asjn[k(r  -  ra)] 
r,  l4-c4AJsin{2*(r-rj))’ 


ra  <  r  <  ra  4- 


In  this  derivation  it  is  supposed  that  the  distortion  is  controlled  solely  by  the  initial  profile;  but  as  propagation 
proceeded,  the  peak  flux  density  would  change  and  the  velocity  excess  would  alter.  There  would  be  a 
tendency  to  increased  flux  density  due  to  bunching  but  this  might  be  outweighed  by  an  opposing  tendency 
represented  by  the  factor  ri/ra  associated  with  spherical  divergence.  Since  the  total  deformation  called  for 
is  only  about  1.5  years  or  so  in  22,  it  has  not  seemed  worthwhile  in  a  first  investigation  to  refine  the  theory 
for  the  waveshape.  By  taking  excess  velocity  to  be  proportional  to  buoyancy  forces,  as  though  viscosity 
were  the  controlling  resistance  to  motion,  we  reach  results  that  will  be  seen  below  to  be  satisfactory.  The 
proportionality  constant  c  will  in  time  be  derivable  from  theory,  pending  development  in  assessing  the 
magnitudes  of  ri  and  BmM(ri)  and  conversely  will  constrain  the  range  of  these  magnitudes.  Meanwhile,  c 
has  been  adjusted  to  meet  a  coarse  requirement  of  1.5  years  advance  in  the  vicinity  of  100  to  150  for  Rm** 

The  general  effect  of  buoyancy  in  a  magnetised  plasma  somewhat  resembles  the  waveshape  effect  on 
intense  sound  waves  or  shallow-water  waves,  where  the  crests  travel  faster  than  the  nodes  and  produce 
wavefront  steepening.  However,  unlike  the  troughs  of  the  mechanical  waves,  which  are  retarded,  both  in- 
phase  and  antiphase  semicycles  of  a  magnetic  wave  are  accelerated.  The  calculated  profiles  of  Fig.  8  show 
the  nature  of  the  phenomenon  described  by  the  simple  buoyancy  theory  presented.  Noteworthy  features  are 
the  wavefront  steepening,  peak  sharpening  and  strengthening,  and  a  break  where  the  descent  from  the  peak 
moderates  into  a  gentler  slope.  For  this  calculation,  k  =  l,  c  =  0.000048,  rx/r3  =  1,  and  A  =  25,  50,  75,  100. 

•  Fig.  6.  A  sinusoidal  profile  A  sin  r  becomes  Asinr/(1  -  cA7  sin2r)  when  buoyancy  is  allowed  for. 
As  the  initial  amplitude  A  passes  through  increasing  values  (25,  50,  75,  100)  the  transmitted  wave  steepens 
in  front,  develops  a  break  on  the  descending  side,  and  the  peak  sharpens  and  intensifies  substantially. 

As  may  be  seen,  the  calculated  profiles  do  not  closely  resemble  actual  sunspot  semicycles  but  this  is 
because  the  three-halves  law  has  not  yet  been  applied.  Fig.  7  (left)  incorporates  this  law  and  also  adjusts  for 
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the  observed  relation  At  =  — 0.042(ilmax  —  72)  for  the  departure  At  in  semicycle  duration  from  the  mean  in 
accordance  with  the  strength  of  the  semicycle  (Bracewell,  1985).  In  addition,  the  maxima  have  been  aligned 
in  time,  for  comparison  with  the  right  hand  part  of  the  figure  (due  to  Waldmeier,  1935). 

The  range  of  semicycle  shapes  provided  by  applying  the  nonlinear  law  and  the  buoyancy  effect  to  semis¬ 
inusoids  is  compared  with  the  full  range  of  historic  semicycles  in  Fig.  8.  The  semicycles  are  arranged  in 
horizontal  rows  according  to  amplitude.  Although  there  is  a  correlation  between  duration  and  amplitude, 
considerable  variation  in  duration  is  apparent  within  each  amplitude  group.  This  is  because  of  other  sig¬ 
nificant  influences  on  duration,  for  example,  strong  odd-numbered  semicycles  in  recent  decades  have  been 
extended  in  duration  by  a  noticeable  negative  offset  in  the  22-year  mean  of  R±(t)-  This  offset  is  analogous 


to  the  additive  350-year  undulation  (Williams,  1981, 1985)  evidenced  by  zig-zag  effect  in  the  rock  record.  So 
a  semisinusoid  was  chosen  to  pass  through  the  start  and  finish  points  of  each  semicycle  and  then  adjusted 
in  accordance  with  Rm*x  for  that  semicycle  by  applying  the  nonlinear  law  to  the  ordinate  and  the  buoyancy 
shift  to  the  abscissa. 

In  addition  to  the  wavefront  steepening  and  peak  sharpening  noted  in  the  stronger  semicycles,  an  initial 
fillet  is  produced  in  semicycles  of  all  strengths.  This  feature  used  to  present  a  difficulty  for  advocates  of  the 
alternating  cycle  because  the  crossovers  were  kinked,  or  inflected  (see  Bracewell,  1953  or  Fig.  2  of  Bracewell, 
1985),  but  now  the  initial  fillet  is  explained  by  the  three-halves  nonlinear  law.  The  break  in  the  descent 
is  seen  to  result  from  buoyancy  in  strong  semicycles  (a  bonus  of  the  buoyancy  calculations,  which  were 
introduced  for  another  purpose).  It  is  noticeable  that  many  of  the  observed  semicycles  have  stronger  breaks 
than  the  simulations,  which  suggests  a  test  for  refinements  to  the  buoyancy  theory. 

It  is  clear  that  the  gross  features  of  the  semicycle  shapes  are  in  accord  with  the  one-parameter  theoretical 
shapes  in  the  majority  of  cases.  The  discrepancies  may  originate  in  various  ways.  First  of  all  the  sunspot 
observations  have  their  own  systematic  but  unknown  errors.  Occasional  semicycles  have  secondary  peaks 
(more  apparent  in  the  monthly  than  in  the  annual  numbers)  which  are  not  simulable  with  one  adjustable  pa¬ 
rameter  Rmax ■  Also,  any  modern-day  350- year  additive  undulation  U(t),  which  would  increase  the  semicvcle 


length  and  enhance  Rmax  where  U(t)  and  R±(t)  had  the  same  sign  (Bracewell  and  Williams,  1986)  has  not 


been  compensated,  nor  has  account  been  taken  of  the  tendency  of  shar 


attributable  to  undulation).  On  the  whole,  the  calculated  shapes,  including  the  year  of  maximum,  agree 
with  the  record  rather  well  within  the  observational  precision  implied  by  the  level  of  irregularities  (“noise”). 

•  Fig.  8.  Comparison  between  recorded  annual  mean  sunspot  numbers  (thin)  and  artificial  profiles 
(thick),  whose  shapes,  including  the  timing  of  the  maximum,  are  determined  by  the  peak  amplitude. 

Eighteenth  century  data 

It  would  take  a  bold  critic  to  claim  that  the  18th  century  profiles,  in  the  left  hand  column,  are  less 
well  represented  than  the  more  recent  profiles.  In  fact,  both  aspects  of  kurtosis,  the  peakiness  of  strong 
semicycles  and  the  inflexion  effect,  as  well  as  the  buoyancy  effect,  are  all  clearly  exhibited,  and  speak  in 
support  of  the  general  reliability  of  the  Zurich  compilation,  confirming  and  extending  an  earlier  discussion 
by  Gleissberg  (1960). 

The  aemicycles  of  1766,  1775,  and  1798  are  reproduced  rather  well.  The  semicycle  starting  in  1755  is 
generally  symmetrical,  but  had  a  peak  in  1761  that  stood  well  above  the  preceding  and  following  years;  this 
peak,  to  which  the  simulated  profile  is  matched,  forces  an  early  simulated  maximum.  If  the  data  reports  are 
erroneous,  one  possibility  is  that  the  sunspot  numbers  for  the  one  year  1761  were  less,  by  about  40  per  cent, 
than  reported.  If  that  single  year  were  altered,  the  simulated  profile  would  fit  the  general  run  well.  However, 
it  is  not  unreasonable  to  accept  that  1761  was  a  spotty  year;  other  years  have  been  seen  when  the  sunspot 
number  exceeded  the  mean  of  the  preceding  and  following  years  by  comparable  amounts.  A  second  possible 
explanation  in  terms  of  erroneous  reports  is  that  the  two  years  1758-1760  were  spottier  than  reported  by 
about  40  per  cent. 

The  .semicycle  starting  in  1784  is  unusual  because  of  its  association  with  the  great  phase  anomaly 
centered  approximately  on  that  year  (Dicke,  1970.  1978;  Bracewell,  1985).  The  reported  maximum  came 
conspicuously  early  but  it  is  hard  to  see  how  this  could  be  due  to  observational  error.  Possibly  the  unveiling 
of  each  semicycle  proceeds  at  a  rate  that  is  advanced  or  retarded  on  occasion  by  convective  motion  in  the 
subphotospheric  medium.  The  shape  described  by  the  reports  is  quite  classical  and  would  closely  fit  a  profile 
lasting  12  years  instead  of  14  years.  The  crossover  year  between  semicycles  4  and  5  is  subject  to  interpretation 
and  in  particular  would  be  subject  to  correction  for  the  350-year  undulation.  The  sense  of  the  correction 
would  be  to  make  the  crossover  year  of  start  later  and  of  finish  earlier.  Perhaps  the  discrepant  simulation  for 
semicycle  4  is  evidence  for  the  presence  of  a  positive  undulation,  which  is  expected,  but  not  clearly  evident 
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(Fig.  2  of  Bracewell,  1985)  in  those  years. 

Origin  of  the  nonlinearity 

Because  the  nonlinearity  discussed  here  was  found  in  both  the  annual  mean  sunspot  series  (Bracewell, 
1985)  and  also  in  the  varve  series  (Bracewell  and  Williams,  1986)  it  is  reasonable  to  look  for  a  physical 
explanation  in  the  sun.  Terrestrial  atmospheric  explanations  could  be  devised  for  either  case  taken  alone  but 
such  explanations  would  call  for  an  influence  of  astronomical  seeing  on  the  counting  of  sunspots  and  sunspot 
groups  on  the  one  hand  and  a  compatible  influence  of  climate  on  varve  thickness  on  the  other.  Turning 
therefore  to  the  sun  we  look  first  at  the  definition  of  the  annual  sunspot  number.  It  is  the  mean  of  the  daily 
values  of  k(lOg  +  /),  where  g  is  the  number  of  groups  and  /  is  the  number  of  spots  counted  by  an  observer 
whose  personal  factor  k  is  chosen  to  establish  historical  consistency. 

There  is  no  reason  to  believe  that  sunspot  number  as  conventionally  defined  should  be  proportional  to 
its  cause  as  measured  in  fundamental  units  of  physics.  It  is  of  course  conceivable  that,  if  the  underlying 
influence  that  causes  sunspots  were  to  double,  then  the  number  of  groups  would  double  in  response.  On  the 
other  hand,  the  groups  might  enlarge  but  not  increase  in  number;  and  since  it  is  known  that  the  number 
of  spots  per  group  increases  rapidly  with  area  of  the  group  (Kiepenheuer,  1953),  here  is  a  possible  source 
of  nonlinearity  built  into  the  sunspot  number  definition.  This  mechanism  could  be  quantified.  Meanwhile, 
there  is  a  different  nonlinearity  that  appears  to  be  sufficient  to  explain  the  three-halves  law. 

Suppose  that  there  is  an  underlying  physical  influence  that  manifests  itself  as  a  center  of  activity  on  the 
solar  surface  and  that  the  “strength”  of  this  center  is  directly  proportional  to  the  influence  I,  by  definition. 
Since  active  centers  erupt  at  intervals  in  various  solar  latitudes  and  longitudes  and  have  finite  lifetimes,  the 
influence  is  sporadic.  Suppose  that  the  ith  sporadic  episode  is  quantified  by  influence  /,•  so  that  £/,-,  summed 
over  a  year,  could  be  the  physical  strength  measure  of  that  year’s  activity.  Now  consider  how  solar  activity 
is  characterized  by  counting  sunspots. 

Each  active  center  is  accompanied  by  a  sunspot  group  whose  temporal  development  is  measured  by 
iv(t),  the  day-by-day  contribution  of  that  sunspot  group  to  the  daily  sunspot  number  Rj(t)  =  £,r,(t). 

Call  a  graph  of  r,(t)  versus  t  in  days  a  sunspot  development  curve.  Such  curves  have  been  much  studied; 
for  example  the  mean  duration  is  6  days  and  (r«(f))  =  12.  Detailed  shapes  of  the  curves  have  been  sorted 
into  categories,  the  sunspot  groups  have  been  classed,  and  the  time  dependence  on  class  has  been  examined 
(  Kiepenheuer,  1953). 

The  annual  mean  sunspot  number  R(t)  is  defined  by  (1/365)  £i*5  Rd(t).  Now  Rj(t)  is  itself  a  sum  over 
the  active  centers;  thus 

36S 

*(<)  =  (1/365)  ££*(*)• 

(=i  • 

By  inverting  the  order  of  summation  we  can  express  this  definition  in  terms  of  the  areas  Ai  =  r,(t)  = 

r<(i),  where  t\,  and  tj,  are  the  day  numbers  for  the  birth  and  death  of  group  i  and  r,(f)  is  deemed  to 
be  zero  before  flt  and  zero  after  t2i.  Thus 

ses 

R(t)  =  (1/365)  ££>(*)  =  (1/365)  £>. 
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The  right-hand  side,  through  the  summation  over  the  groups  »,  represents  a  mean  that  is  weighted,  very 
reasonably,  in  accordance  with  the  frequency  of  each  class  of  sunspot  group;  in  addition  weight  is  given  both 
to  sunspot  size  and  to  longevity  through  the  contribution  to  area  Ai  under  the  development  curve.  We  see 
that  R(t),  while  it  is  certainly  a  measure  of  the  year’s  spottiness,  is  not  at  all  a  measure  of  the  year’s  activity 
as  expressed  by  Lit.  To  take  an  analogy,  £/,  is  like  the  annual  rainfall,  while  £  Ai  is  like  the  annual  mean 
depth  of  rainwater  standing  in  a  pond;  one  heavy  fall  of  rain  will  contribute  more  to  the  mean  depth  of  water 
than  the  same  amount  of  rain  falling  in  showers.  In  other  words,  when  an  episode  is  big  as  quantified  by  /,• 
then  not  only  its  size,  but  also  its  duration  is  factored  into  Ai.  For  this  reason  alone,  R(t)  in  years  of  large 
annual  sunspot  number  will  overestimate  £/<  simply  because,  by  and  large,  big  sunspots  are  more  durable. 

To  estimate  the  quantity  £/<  we  would  need  to  take  all  the  sunspot  group  development  curves  for 
the  year  and  sum  their  amplitudes,  i.e.  their  contribution  to  the  sunspot  number  on  the  day  when  that 
contribution  to  the  sunspot  number  was  a  maximum.  However,  since  much  statistical  analysis  is  on  record, 
we' can  simply  plot  the  areas  A  of  different  sunspot  development  curves  against  Rm^c  the  maximum  sunspot 
number  reached,  where  the  bars  indicate  averages  taken  within  the  different  categories.  The  results  deduced 
from  Kiepenheuer  (1953)  are  plotted  in  Fig.  9. 


•  Fig.  9.  The  number  of  sunspot-days  A  representing  the  mean  area  under  different  categories  of 
sunspot  group  development  curves  as  a  function  of  Rm^x,  the  maximum  contribution  to  sunspot  number 
reached  in  the  course  of  development  averaged  over  the  same  category.  The  straight  line  has  a  slope  of  3/2. 

We  see  that  on  a  log-log  plot  of  A  versus  Rmxx  the  12  points  lie  near  a  straight  line  whose  equation  is 

A  =  ZR*Z- 

This  result  could  be  refined  by  weighting  the  points  in  accordance  with  the  frequency  distribution  of  the 
different  categories  of  development  curves.  Meanwhile,  it  is  extremely  striking  that  the  same  three-halves  law 
emerges  directly  from  the  available  statistical  studies  as  was  arrived  at  from  the  third-harmonic  investigations. 
One  concludes  that  the  gross  phenomenon  of  overcounting  sunspots  in  accordance  with  the  duration  of  their 
groups  is  sufficient  to  explain  the  bulk  of  the  nonlinear  phenomenon  that  was  discovered. 

It  ought  to  be  added  that  the  linearized  variable  Ru„ ,  that  has  been  di&  -*d  by  terms  such  as  “physical” 
and  “underlying”,  and  which  is  the  goal  of  a  deliberate  quest  for  something  “fundamental,”  by  no  means 
damages  the  physical  status  of  the  integrated  quantity  R(t),  which  takes  account  of  duration.  The  latter  is 
after  all  the  quantity  which  correlates  with  ionospheric  characteristics  for  the  simple  reason  that  sustained 
solar  ultraviolet  radiation  from  a  long-lived  centre  of  activity  continues  to  produce  ionization  as  long  as  it 
is  present.  Assume  that  varve  thickness  represents  the  annual  run-off  from  a  periglacial  region  into  a  lake, 
and  that  the  run-off  in  turn  measures  the  integrated  solar  ultraviolet  which,  not  being  impeded  by  today’s 
ycygen  atmosphere,  could  reach  down  to  the  troposphere  in  Precambrian  times.  Then  it  is  the  sunspot 
number  R(t),  not  Riin(t),  with  which  varve  thickness  would  be  expected  to  correlate:  the  presence  of  the 
nonlinear  effect  in  varve  thickness  is  consistent  with  this  expectation.  Similarly,  annual  rainfall,  which  is  like 
R(t),  and  the  nonlinearly  related  mean  depth  of  a  rain-fed  pond,  which  is  like  Run(i),  are  each  valid  in  their 
own  applications. 

Sunspot  area 

Authors  have  often  spoken  respectfully  of  sunspot  area,  a  measurement  that  possesses  long  term  stability, 
while  nevertheless  using  sunspot  number  because  it  has  been  available  in  real  time,  whereas  the  Greenwich 
Phoioheltographic  Results  used  to  appear  many  years  late.  For  the  years  1784  to  1938  it  was  found  that 
sunspot  area  F,  measured  in  millionths  of  the  solar  hemisphere,  was  related  to  the  sunspot  number  R 
by  F  =  16.7 R.  But  in  fact  F  is  not  proportional  to  R.  Data  from  both  Greenwich  and  the  U.S.  Naval 
Observatory  (reproduced  by  Waldmeier,  1968)  respectively  fit  the  following  power  laws:  F  =  1600(/2/120)1J5 
and  3400(72/180)1J  (Fig.  10).  Consequently  sunspot  area  should  be  proportional  to  R}^  approximately 
and  exhibit  even  greater  departure  from  the  underlying  influence  than  does  the  sunspot  number  itself.  The 
explanation  for  the  nonlinear  relation  between  sunspot  area  and  sunspot  number  clearly  lies  with  the  way  in 
which  sunspot  groups  of  different  total  area  are  divided  into  separate  spots.  One  could  rederive  the  linearized 
sunspot  number  Run(t)  from  sunspot  area  records,  and  possibly  obtain  a  better  result  than  working  from 
sunspot  numbers  allows. 

•  Fig.  10.  Data  for  annual  mean  sunspot  area  F  versus  annual  mean  sunspot  number  A  (Greenwich 
below,  Washington  above)  showing  a  clearly  nonlinear  relation. 


Conclusion 

A  nonlinearity  entering  into  the  definition  of  sunspot  number  has  been  demonstrated  and  explained. 
For  purposes  of  physical  theory  it  will  be  desirable  in  future  to  correct  sunspot  numbers  for  the  nonlinear 
relationship  that  is  now  known  to  exist.  As  an  example  of  such  a  correction,  combined  with  allowance  for 
buoyancy  in  a  magnetized  plasma,  it  has  been  possible  to  reproduce  theoretically  the  known  Waldmeier 
relation  between  rise  time  and  maximum  sunspot  number,  and  other  features  of  the  11-year  sunspot  curve. 
A  break  in  the  decline  of  strong  11-year  cycles,  which  has  been  demonstrated  to  exist,  is  predicted  by  the 
same  theory.  Eighteenth  century  sunspot  numbers  fit  theoretical  semicycle  shapes  quite  well;  their  reputed 
unreliability  is  not  confirmed  by  the  present  work.  The  physical  model  envisages  an  upward  travelling  22- 
year  wave  launched  by  a  reciprocating  magnetic  oscillator  deep  in  the  sun.  While  this  model  is  at  present 
incomplete,  it  is  qualitatively  compatible  with  the  fast  rise  times  associated  with  strong  11-year  cycles, 
and  with  other  features  of  the  cycle  shape  discussed  in  earlier  papers,  and  to  this  extent  is  supported  by 
observation. 
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•  Figure  1.  Features  of  the  shape  of  a  strong  semicyde  (left)  indude  a  leading  fillet,  a  convex 
segment  peaking  early,  and  a  relatively  long  decline,  sometimes  with  a  noticeable  break  in  slope. 
A  weak  semicyde  (right)  tends  to  be  much  more  symmetrical. 


•  Figure  2.  A  nonlinear  characteristic  (top  left)  relates  a  stimulus  x(t)  (lower  left)  to  a  re¬ 
sponse  y(t)  (upper  right).  In  terms  of  the  corresponding  points  (xm,  tm)  and  (ym,  tm)  the  reduction 
factor  K  is  xm/ym. 


•  Figure  3.  The  nonlinear  relation  between  R  dk  (!)  and  Riin(t)  pot  in  evidence  by  yearly 
point*  from  1700  to  1985.  The  curve  is  R±(t)  =  100[|  fin„(t)  |  /83]3/,agn.RJi„. 


•  Figure  4.  The  data  of  the  previous  figure  combined  in  one  quadrant  by  plotting  log  |y| 
against  log  |*|.  The  straight  line  has  a  slope  3/2. 


#  Figure  5.  Artificial  “11-year*  semicyd*  shape*  versus  time  (top)  associated  with  fixed- 
duration  semisinuoids  (center)  of  four  different  strengths.  The  artificial  shapes  are  derived  from 
a  three-halves  nonlioear  law  without  introduction  of  asymmetry.  The  difference  curves  (below) 
resemble  third  harmonic,  in  the  middle  range  of  strength. 
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•  Figure  6-  A  sinusoidal  profile  Asinr  becomes  Asinr/(1  —  cA2sin2r)  when  buoyancy  is 
allowed  for.  As  the  initial  amplitude  A  passes  through  increasing  values  (25,  50,  75,  100)  the 
transmitted  wave  deepens  in  front,  develops  a  break  on  the  descending  side,  and  the  peak  sharpens 
and  intensifies  substantially. 
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abscissae :  years  from  maximum 
ordinates:  relative  numbers 


•  Figure  7.  Artificial  semicycle  shapes  incorporating  (a)  nonlinearity,  (b)  an  advance  in  max¬ 
imum  associated  with  magnetic  buoyancy,  and  (c)  a  shortening  of  the  semicycle  duration  inversely 
correlated  with  semicycle  strength  (left).  For  comparison,  some  selected  semicycles  from  history 
(right).  Note  that  the  historical  semicycles  fail  to  reach  zero  because  spots  from  the  preceding  or 
following  semicycles  are  included. 
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•  Figore  8.  Comparison  between  recorded  UMtl  mu  soaspot  ttmbwi  (thin)  sad  artificial 
profiles  (thick),  whose  shapes,  iadadiag  the  timing  of  the  maximum,  sis  determined  bp  tbs  peek 
amplitude. 


e  Fig* re  9.  The  somber  of  sunspot-days  A  representing  the  mesa  use  onder  different 
categories  of  soaspot  group  development  corves  a*  a  function  of  La,  the  maximum  contribotioe 
to  soaspot  aomber  reached  ia  the  coon*  of  development  averaged  over  the  same  category.  The 
straight  line  hat  a  slope  of  3/2. 


